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The relation between the long wavelength limit of solutions to the cosmological perturba- 

^ ' tion equations and the perturbations of solutions to the exactly homogeneous background 

■*>>> ■ equations is investigated for scalar perturbations on spatially flat cosmological models. It 

Q-i| is shown that a homogeneous perturbation coincides with the long wavelength limit of 

^ ', some inhomogeneous perturbation only when the former satisfies an additional condition 

^' corresponding to the momentum constraint if the matter consists only of scalar fields. In 

contrast, no such constraint appears if the fundamental variables describing the matter con- 

^ ■ tain a vector field as in the case of a fluid. Further, as a byproduct of this general analysis, 

it is shown that there exist two universal exact solutions to the perturbation equations in 

the long wavelength limit, which are expressed only in terms of the background quantities. 

They represent adiabatic growing and decaying modes, and correspond to the well-known 

exact solutions for perfect fluid systems and scalar field systems. 
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§1 Introduction 

In the current standard scenario based on the gravitational instabihty theory, the present 
large scale structures of the universe are formed through the following four stages: produc- 
tion of seed fluctuations in the early universe, their linear evolution on superhorizon scales, 
the subsequent linear modulation after they enter the horizon, and the final non-linear 
evolutions. In the inflationary universe models the seed fluctuations in the first stage are 
produced from quantum fluctuations on the Hubble horizon scales during inflation, and 
we now have universal formulae to determine their nature such as the amplitudes and the 
spectrum for a wide variety of inflation models. We can now also easily handle the evolu- 
tion during the third stage because the matter content of the universe during this stage is 
rather simple and restricted. Of course the actual behavior of perturbations on subhorizon 
scales are quite complicated even in the linear stage, and we need numerical computations 
to determine their details. 

On the other hand universal formulae to determine the evolution of perturbations during 
the second stage have not been established yet, although it is generally believed that the 
so-called Bardeen parameter is conserved with good accuracy during this stage and this 
conservation law essentially determines the amplitudes and the spectrum of perturbations 
when they reenter the horizon, from those at the first stage ||l|]. Of course the conservation 
of the Bardeen parameter during the Friedmann stage is well established including the case 
in which the equation of state of cosmic matter changes slowly 0, ^. The simplicity comes 
from that fact that the evolution of perturbations on superhorizon scales are determined 
with good accuracy by that on the long wavelength limit 0, |^, ^ . In this limit the evolution 
of adiabatic modes is determined by a set of simple Ist-order ordinary differential equations 
for two gauge-invariant variables, from which follows the exact conservation of the Bardeen 
parameter |^, ^. However, for some important situations, such as the reheating stage 
of inflation models, this simple analysis does not work because the adiabatic mode may 
produce entropy modes which feed back to the behavior of the adiabatic modes[0. In 
these situations the knowledge on the evolution of perturbations of all the components of 
matter in the long wavelength limit is required to determine the behavior of the Bardeen 
parameter. 

In this connection Nambu and Taruya recently wrote an interesting paper p] in which 
they stated that solutions to the gauge- invariant perturbation equations for a multi-component 
scalar field system on expanding universe in the long wavelength limit are obtained as 
derivatives of exactly homogeneous solutions to the Einstein equations with respect to the 
solution parameters. If this result is correct and all the solutions to the perturbation equa- 
tions can be obtained in this way, it provides a very powerful tool to analyze the behavior 
of perturbations during inflation and reheating including the problem of the conservation 
of the Bardeen parameter. However, no proof nor explanation on their statement was given 
in their paper. Further, direct calculation shows that the perturbation equations in the 
long wavelength limit do not coincide with those for exactly homogeneous perturbations. 

In order to see whether Nambu and Taruya's statement is true or not as well as to 
obtain a deeper understanding on the behavior of perturbations on superhorizon scale per- 
turbations, in the present paper, we investigate the relation between the solutions to the 
perturbation equations in long wavelength limit and the exactly homogeneous solutions 
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to the Einstein equations for multi-component systems on spatially flat Robertson- Walker 
universe. Our arguments are quite general except for the assumptions that the fundamental 
variables describing matter are scalar fields and/or vector fields, and that the amplitude of 
anisotropic stress perturbations vanishes rapidly enough in the long wavelength limit. We 
will clarify general conditions under which the perturbation solutions in the long wavelength 
limit coincide with some exactly homogeneous perturbations to the homogeneous solutions 
to the Einstein equations. In particular, we will show that for the mult i- component scalar 
field system not all the exactly homogeneous perturbations directly correspond to pertur- 
bation solutions in the long wavelength limit. Further, as a byproduct of our analysis, we 
will show that the perturbation equations in the long wavelength limit have two universal 
adiabatic solutions which can be expressed explicitly as time integrals of known background 
quantities. 

The paper is organized as follows. First in the next section we give basic definitions of 
perturbation variables used in the present paper and the perturbation equations for them. 
By inspecting the dependence of these quantities and equations on the wave number k of 
perturbations, in §3 we derive the conditions on the exactly homogeneous perturbations 
to coincide with the fc — > limit of some solutions to the perturbation equations with 
k y^ 0. In §4, as a special application of the argument in §3, we show the existence of 
two universal adiabatic modes in the k -^ limit. Then in §5 and §6 we specify our 
arguments to a multi-component scalar field system and a multi-component perfect fluid 
system, respectively, and clarify the relation between the perturbations in the k —>■ limit 
and those with k = 0. Section 7 is devoted to summary and discussions. 

§2 Perturbation Equations 

In this section we recapitulate the definitions of basic perturbation variables and their 
equations in the framework of the gauge-invariant perturbation theory. We adopt the 
notations used in the review article by Kodama and Sasaki 0. 

We only consider perturbations on a spatially flat(i^ = 0) Robertson- Walker universe 
throughout the paper. Hence the background metric is given by 

ds^ = -dt^ + a{tfdx^, (2.1) 

and its perturbation by 

ds'^ = -(1 + 2AY)dt^ - 2aBYjdtdx^ + a^ [(1 + 2HLY)6jk + ^HrYjk] dx^dx\ (2.2) 

where Y, Yj and Yjk are harmonic scalar, vector and tensor for a scalar perturbation with 
wave vector k on fiat 3-space: 

Y:=e^'-, Yr-=-^^iY, F^., := (^1^^., - M^) K (2.3) 

Under the infinitesimal gauge transformation 

St = TY, 6x^ = LY\ (2.4) 



the metric perturbation variables in Eq. (|2.2|) transform as 

A = A-f, (2.5) 

k 
B = B + aL + -T, (2.6) 

a 

Hl = Hl-^L~ HT, (2.7) 

Ht = Ht + kL (2.8) 

where the overdot represents the derivative with respect to the proper time t and H is the 
cosmic expansion rate a/a. From this it follows that the spatial curvature perturbation TZ 
and the shear ag defined by 



n:=HL + ^Ht, 


(2.9) 


ag:=^HT-B 


(2.10) 


transform as 

n = lZ-HT, a„ = a„- -T. 

a 




(2.11) 


Hence we obtain the following two independent gauge-invariant combinations: 




A:=A-{n/Hy, 


(2.12) 




(2.13) 



All the other gauge-invariant combinations of metric perturbation variables are written as 
linear combinations of A and $. 

Here note that for the exactly homogeneous perturbations corresponding to fc = 0, the 
vector-like and the tensor-like perturbation variables should vanish. Hence for the metric 
perturbation, B = Ht = 0. Correspondingly, the gauge freedom for them is just the 
time reparametrization T, and L should vanish in the gauge transformation. Therefore 
there exists only one gauge-invariant combination for the exactly homogeneous metric 
perturbations, which coincides with A. The variable $ has no counterpart for them. This 
point is very important in the argument on the correspondence between the perturbations 
in the k —>■ limit and the exactly homogeneous ones. 

Next we consider the matter perturbations. In this paper we assume that the funda- 
mental variables describing matter are scalar and vector quantities. This assumption is 
satisfied in most of the realistic applications. Let us denote the corresponding perturba- 
tion variables by xi ^ind vp where I and P are indices labeling components. Then they 
transform under the gauge transformation as 

Xi = Xi - SiT, vp = vp + at, (2.14) 

where Sj is the background quantity corresponding to the perturbation xi- Hence the 
corresponding gauge-invariant variables are given by 

Xi := XI - §7^, (2.15) 

Vp := Vp - ^Ht. (2.16) 

k 



In the exactly homogeneous system the matter is described only by the scalar quantities 
Sj. Hence its perturbation gives only the scalar-type gauge-invariant variables Xj as in 
the case of the metric perturbation. 

The Einstein equations relate these matter variables with the metric variables through 
the energy-momentum tensor. For scalar perturbations its generic form is written as 



rfiO 

rpO 

j 
fj 



a{p + p){v-B)Y,, 
' ^6p6i + UY,' 



k li 



(2.17) 
(2.18) 
(2.19) 



where p and p are the background values of the energy density and the pressure, respec- 
tively, and follow the equations, 

,^2 



E' 



-P, 



p = -3(p + p)i7. 



(2.20) 
(2.21) 



In the present paper we only consider the case in which the anisotropic stress perturbation 
n, which is gauge-invariant by itself, vanishes faster than k^ in the fc — ;> limit. 

Applying the general argument on the matter perturbation above, we can construct 
from the density perturbation (5p, the velocity perturbation v and the isotropic stress 
perturbation 5p the following three gauge-invariant combinations: 



pA3:=5p + 3(p + p)7^, 

V ■.= v 






pT := Sp — c^Sp 



(2.22) 
(2.23) 
(2.24) 



where cj. := p/p is the square of the sound velocity. 

Now in terms of these gauge-invariant variables the perturbation of the Einstein equa- 
tions are written as 






^'T,^ ■■ 



A 



GO 






Gk 



G\S- 



k 



.-Tl. 



n-Tj: 



^\n 



k 



2 ^ 
3 



3a2F2 



A+-il + w)Z 



0, 



T5' 



«xi^ 



I'-'k) 



HA + {2H + 3H^)A 



K 



{pT + clpA, 



K 



A+-(-^ 
a \H 
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where w = p/p and Z is the Bardeen parameter defined by 

aH „ aH , 

Z := $ - ^-V = 7^ - ^{v - B) 



k k ' 

If we eliminate A^ from Eq.( |2.27l ) using Eq.( |2.25| ), we obtain 



A 



1 



w 



1 e 



l + wHa? 
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$ + 



H 



p + p V2 



-pV -li 



(2.25) 
(2.26) 

nj(2.27) 
(2.28) 



(2.29) 



(2.30) 



Note that we have not multiphed any power of k in deriving these equations from 
the Einstein equations. Hence putting k = in these equations yields the equations for 
the exactly homogeneous perturbations if one takes account of the fact that the terms 
multiplied by inverse powers of k vanishes identically for these perturbations as mentioned 
above. 

In general these equations must be supplemented by the expressions for A^, Z, T and 11 
in terms of the fundamental gauge-invariant matter variables Xj and Vp and their evolution 
equations. We do not write them explicitly here because the details of these equations are 
not necessary until we specialize the general arguments to specific models. 

§3 /c = vs /c ^ limit 

In this section we clarify under what conditions the k ^ limit of solutions to Eqs.( p^.25D - 
(|2.28| ) coincide with some solutions to the corresponding equations for k = 0. For that 



purpose first note that the gauge-invariant variables introduced in the previous section are 
classified into two groups. 

The first one consists of the gauge-invariant variables whose expressions in terms of the 
gauge-variant perturbation variables do not contain inverse powers of k. A, Ag, T and Xj 
belong to this group. These variables always have well-defined finite fc — > limits which 
coincide with the corresponding gauge-invariant variables for the exactly homogeneous 
perturbations. Of course in taking k ^ limit the conditions 

B,Ht,vp,v,U/P ^0 (3.1) 

should be taken into account. 

On the other hand the second group consists of the gauge-invariant variable which 
contain inverse powers of k when expressed in terms of gauge-variant variables. $, Z{ot 
V) and Vp belong to this group. These variables have no counterpart for the exactly 
homogeneous perturbations, and their A; ^ limits depend on how fast the quantities in 
Eq.( |3.1|) vanish in the limit. 

With these points in mind, let us first consider the fc — ^ limit of the perturbation of the 
equations for the fundamental matter variables. When expressed in terms of the original 
gauge-variant perturbation variables x/, vp, A, B, H^ and Ht, they do not contain inverse 
powers of k. We can write them in terms of gauge-invariant variables in the following way. 
First eliminate xi ^^^^ ^ by replacing them by Xj, A and TZ. Second eliminate vp using 
VpQ = vp—VQ and Z. These procedures do not produce terms with inverse powers of k, and 
the coefficient of Z contains a positive power of k. Now since all the gauge-variant matter 
variables and A are eliminated, the remaining terms should be proportional to $( and its 
time derivatives). This implies that in the final expressions Z and $ appears only in the 
forms k"^Z{m > 1) and A;"<l>(n > 2). Further the other terms, which are proportional to 
A, Xj or VpQ, contain no inverse power of k. Therefore, taking into account the condition 
(|3.1|), the fc — *■ limit of a solution to the perturbation equations satisfies the corresponding 



equations for /c = if and only if 

fc2$ -^ 0, (3.2) 



kZ -^ 0, (3.3) 

VpQ -^ 0. (3.4) 

Next let us examine the Einstein equations. Among the four equations it is easy to see 
that Eq.(|]25D and Eq. (g^) reduce to the perturbation of Eq.( ^^ and Eq.(|]2l]) in the 



/c — !► hmit under the condition fc^$ — >■ 0. On the other hand Eq.( |2.26D and Eq.( |2.28D 
become trivial for the exactly homogeneous perturbations. However, for /c 7^ 0, these 
equations yield non-trivial perturbation equations 

^+^(l + w)Z = 0, (3.5) 

1 / a 



A^-y-^) =0, (3.6) 

where we have put 11 = because it does not affect the arguments on the fc ^ limit 
under the assumption II/A;^ -^ adopted in the present paper. 
The second equation of these can be solved in terms of $ as 

^ = — (c-f a{t)A{t)dt] , (3.7) 

where C is an integration constant and to is an initial time. This equation can be regarded 
as an equation to determine the k ^ limit of $ from a solution for A to the perturbation 
equations with k = 0. In this viewpoint the consistency condition (|3.2| ) is simply replaced 
by the condition on the fc-dependence of C, 

k^C{k) -^ 0. (3.8) 

On the other hand the condition (|3.3|) is always satisfied under Eq. (p.5|) . Further the 
condition ( |3^ ) just select a subset of solutions to the perturbation equations with A; 7^ 0, 
and does not give any restriction on the exactly homogeneous perturbations. Therefore 
starting from any solution to the perturbation equations with k = 0, one can always 
construct the gauge-invariant quantities representing the k —>■ limit of a solution to the 
perturbation equations with /c 7^ by supplementing A and Xj for the exactly homogeneous 
perturbation with $ determined by Eq. (|3.7| ), provided that Eq.( p.5D do not yield any 
additional constraint on the seed exactly homogeneous perturbation. 

In the case in which the fundamental variables describing matter contains a dynamical 
vector field, the k ^ limit of Z depends on the value liTak-,o{v — B)/k which cannot 
be determined from the information of the exactly homogeneous perturbations. Hence 
Eq.( p.5|) can be simply regarded as the equation to determine the k ^ Q limit of Z. 

In contrast, in the case in which the matter is described only by scalar fields, v — B 
in T? should be written as a combination of the spatial derivatives of xi- Hence it must 
be proportional to a positive power of fc, and Z has a well-defined /c — *■ limit which is 
written only in terms of Xj. Thus the condition Eq.( p.5|) yields a restriction on the seed 
exactly homogeneous perturbation in order for it to be a /c —^ limit of some solution to 
the perturbation equations with fc 7^ 0. 



14 Universal Adiabatic Solutions 



Let a(t) and Si{t) be an exactly homogeneous solution describing the background. Then, 
since the scale factor a{t) conies into the evolution equations only through H = a/a, 
(1 + X)a(t) and Si(t) also gives an exactly homogeneous solution to the Einstein equation 
where A is a constant. If we regard this solution as an exactly homogeneous perturbation, 
all the corresponding gauge- variant perturbation variables for matter vanish, and the metric 
perturbation is given by 

A = B = 0, n = HL = \. (4.1) 



From this we obtain 



A 



[l + w)\. 



Equation ( p.7| ) determines $ from this as 



$ = C— + ^A— / (1 + w)a{t)dt. 
a 2 a Jto 



(4.2) 



(4.3) 



If there exists a dynamical vector field describing matter, this equation and Eq. 
determines Z as 

Z = —^ — r^=A. (4.4) 



3(1 



w] 



On the other hand, if the matter is described only by scalar fields, v — B is written as 
V — B = kx where x is a combination of the scalar field perturbations xi which does not 
contain a negative power of k. Since the matter perturbation vanishes in the present case, 
this term should vanish, which implies that Z = X — aHx = A. Hence the condition ( p.5| ) 
is satisfied. 

Thus we find that the perturbation equations have always two solution for which the 
fc — > limits of A, $ and Z are given by Eqs. (|4.2|) - ([4.4| ). Since the matter is not perturbed, 
r = for them. Hence they represent adiabatic modes in the fc ^ limit. Clearly the solu- 
tion proportional to C is a decaying mode, while that proportional to A is a growing mode 
because the Bardeen parameter Z is a non- vanishing constant. Note that these solutions 
are universal in the sense that they are valid for any matter contents and interactions. 

This universal solutions cover most of the exact solutions in the A; ^ limit found so 
far, and the general argument developed here explains why such exact solutions were found 
for large variety of matter content. For example, for the case in which matter consists of a 
single-component perfect fluid, there exist two adiabatic modes. The above solutions just 
give their /c ^ limit, and in terms of the standard variables A = 6p/ p + 3(1 + w)aH{v — 
B)/k and V they are expressed as 



A 
V 



2k' 



.$ 



2P 



ah a^ 



3 f 

C + -\ {l + w)adt 
2 Jto 

C-X (^a{to)/H{to) + 



adt 



to 



(4.5) 
(4.6) 



which recover the well-known exact solutions! 



Another example is the case in which matter consists of a muhi-component scalar field 
= (0/). In terms of the gauge-invariant variable for the scalar field perturbation defined 
by 

Xi = 5<Pi - ^7^, (4.7) 

the above solution gives 

Xi = -^A, (4.8) 

and 

i = C^^4^[p^. (4.9) 

a 2 a J H^ 

These give the extension of the well-known exact solution in the fc — > limit for the single 
component scalar field casef^, ||, |^. 



1 5 Scalar Field Systems 



In this section we examine how the condition ( p.5| ) is expressed explicitly for a multi- 
component scalar field system. 

The Lagrangian density for a multi-component scalar = (0/) is generally expressed 

as 

^g'-'d,^ ■ 9.0 + f/(0) 



C = -V^ 



(5.1) 



where f/(0) is a potential. For a homogeneous background the energy density p and the 
pressure p are expressed as 

p=^02 + f/(0), p=l02_f/(0)^ (5.2) 

where 0^ = ■ 0. The background equation of motion and the Einstein equations are given 
by 

^ + 3H(j) + DU = 0, (5.3) 

where DU = {dU/d(j)i), and Eq.dOq). 

Since the energy-momentum tensor for the Lagrangian density (|5.1| ) is given by 



Ti: = V^0 ■ 9,0 - i(5^ (V> ■ Va + 2U) , (5.4) 

the perturbation variables for the energy-momentum tensor are expressed as 

Sp = -A02 + 050 + DU ■ 6(j), (5.5) 

{p + p){v - B) = -<j) ■ 6(j), (5.6) 
a 

6p = -A02 + 0(50 - DU ■ 5(1). (5.7) 



Hence A^ and Z are expressed in terms of the gauge-invariant variable X for the scalar 
field perturbation defined by Eq. (|4.7|) as 

(5.8) 

(5.9) 



pA, = -A(f + (j)-X + DU-X, 



-H- 



■X 



Inserting this expression for A^ into Eg. ( p. 25 ), we obtain 

1 fc2 



2UA + (j)-X + DU ■X = 2^'—^. (5.10) 

From this, the expression for Z and the background field equation, it follows that 

^,5(l,„,i^.-^.W|.A-).^^*, (5.11) 



where 

iy(Xi,X2):=Xi-X2-Xi-X2. 

Hence the condition ( |3.5|) is expressed as 



(5.12) 



limW^ — ,X =0. 



(5.13) 



Now let us see what kind of restriction this condition gives for exactly homogeneous 
perturbations. Let us denote the gauge-invariant variable for an exactly homogeneous 
perturbation 6(f) of the scalar field a.s X = 5(j) — TZ4>/H, where 



TZ = Hl = 5a/ a. 
Then from the perturbation of Eq.( |5.3|) and Eq. (|2.20|) , we obtain 



X + ?,HX + D'^U{X) ~4)A + 2DUA = 0, 
2UA + d)-X + DU-X = 0, 



(5.14) 



(5.15) 
(5.16) 



where D^U{X) = {Xjd^U/d(f)id(j)j). Here note that X = ^/H and A = 30V(3p) = 
3(1 +w)/2 is a solution to this equation, which corresponds to the adiabatic growing mode 
obtained in the previous section. 

Eliminating A from these equations we obtain the following second-order differential 
equation for X: 

LiX) = -!^(l]w(lA, (5.17) 



U \H. 



W 



where 



L{X) ■.= X + 3HX + 



H? ( c? 



D'V--A-^.. 



X. 



(5.18) 



From this it follows that W = W{(j)/H, X) follows the differential equation 

W -UwH + yDU -^W = Q. (5.19) 

This equation shows that if W vanishes at an initial time t = t^, it vanishes at any 
time. Hence Eq. (|5.13| ) reduces to the condition on the initial value of solutions to Eq.( ^.17l ). 



Therefore, taking into account the fact that X = (p/H is one of such solutions, we find 
that for a A^-component scalar field system {2N — 1) independent solutions to Eq.( ^.17[) 
correspond to the /c — *> limit of solutions to the perturbation equation with fc 7^ 0. 

Here note that we have already obtained two universal solutions in the k ^ Q limit in 
the previous section. Among them, the growing mode, coincides with one of the solutions 
obtained from X . On the other hand the decaying mode is not contained in the latter 
because X = {5a/ a = 0) for it. Hence we have obtained 2N independent solutions in the 
/c — > limit. Since the number of the dynamical degrees of freedom of the A^-component 
scalar field system is 2N, these exhaust all the solutions to the perturbation equations in 
the k —>■ limit. However, if one wants to determine the fc-dependence of the solutions 
around k = beyond the k —>■ limit by solving the perturbation equations iteratively 
with respect to k, one needs all the solutions to the perturbation equation for X in the 
/c — i> limit to find the Green function. Fortunately we can obtain the remaining one 
independent solution for X explicitly from a solution for X with W ^ in the following 
wayQ. 

First note that the perturbation equation for X with A; 7^ is expressed in terms of the 
operator L defined in Eq.( p.l8D as 

1.2 
L{X) + — X = 0. (5.20) 

From this we see that in the A; — > limit X satisfies the equation L{X) = 0, which confirms 
the above argument. 

On the other hand from Eq.( |5.19| ) it follows that 



From this find that 



^-^W (^,x]= const. (5.21) 



^^^'^/^^^i^^^'^^'+^Ai ^^-^^^ 



coincides with the right-hand side of Eq.( |5.17| ) if / satisfies 
This implies that 



•^This point is first suggested out by Sasaki and Tanaka from the consideration on the role of the /c — > 
hniit of Ht in the discrepancy between the equations for X and A" (private communication) 
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satisfies L{X) = for any solution X to Eq. (|5.17|) . Since tlie second term on tlie riglit- 
liand side of tliis equation vanishes if W{(j)/H, X) = 0, and W{(p/H, X) does not vanishes 
iiW{^/H,X) ^Ofrom 

X given by Eq. (|5.24|) exhausts all the solutions to the perturbation equation for X in the 
A; — > limit. 

§6 Perfect Fluid Systems 

In this section we apply the argument in §3 to a multi-component perfect fluid system as 
an example of non-trivial systems in which matter variables contain a djTiamical vector 
field. 

The equations of motion of a perturbed mult i- component system are given by 

V.f,; = Qj^ = Qju^ + fi^, (6.1) 

where Qj^ represents the energy-momentum transfer term for the component /, u'^ is the 
4- velocity of the whole matter system, and Qj := —u^Qifj,- Because of the conservation of 
the total energy-momentum, Qj^ satisfies 

E^/M = 0. (6.2) 

For the scalar perturbation, the energy-momentum tensor and the energy-momentum 
transfer vector of each individual component are expressed as 

fj',=-ipj + 6pjY), (6.3) 

fj!; = a{pj+pj){vi-B)Y^, (6.4) 

fji = {pj6i + 6pi6i + UjYi), (6.5) 

Qio = -[Qi + {QiA + 6Qi)Y], (6.6) 

Qij = a[Qiiv -B) + F,j]Y,, (6.7) 

where pi, pi and Qi are the background values of the energy density, the pressure and 
the energy transfer of the component /, respectively. In terms of these quantities the 
background part of Eq. (|6.1j) is written as 



pi = -3hj{l-qj)H, (6.8) 

where 

hi := PI + PI, qi := Qi/{3Hhi). (6.9) 

As in §2, we construct from the density perturbation 6pi, the velocity perturbation vi, 
the isotropic stress perturbation 6pi and the energy transfer perturbation 6Qi the following 



11 



four gauge-invariant combinations: 



(5p/ + 3(p/+p/)7^, 



Vi ■= vi - -Ht, 
piTj ■= Spi - cjSpi, 



QiE, 



9I 



SQi - ^7^ 



(6.10) 
(6.11) 
(6.12) 

(6.13) 



where cf = pi/ pi- The anisotropic stress perturbation I!/ and the momentum transfer 
perturbation F^i of the component / are gauge-invariant by themselves. 

/^From the relation T^,y = Yl,iTi^y the variables for the whole system are expressed in 
terms of those for each component as 



I I 

pAg = Y^ PlAgJ, hV 

I 



Y^hjVi, pT = Y.piTj + pTrei, n = ^n,, (6.15) 
II I 



where pT 



rel 



Y.iicj — c'^)piAgj. Further Eq.( |6.2| ) gives the constraints 

^Q, = 0, Y.QiEgi = o, E^c/ = o. 



(6.16) 



In terms of these gauge- invariant variables the perturbed equations of motion Eq.( |6.1| ) 
are written as 



Vv-f/o 



Qio 

^uTj^ = Qii : 



^2 

{piAgiY + SHpiAgr + -YTrhii^ - Zi) + ^H{piTi + c]piAgi) 

an 



QiA + QiEgi 
'hiZj\' 



H 



where 



-aFci + k—rZ 
rl 



Zi:=^-^Vi 



ShiZi + hi A + piTi + CjpiAgi - -11/ 



Qi 



n 



fivi-B). 



(6.17) 



(6.18) 



(6.19) 



In deriving these equations from the equations of motion Eq. (|6.1|) , we have not multiplied 
any power of k. Hence putting /c = in these equations yields the equations for the exactly 
homogeneous perturbations. 

Now let us examine the relation of the fc — ;> limit of these equations and the corre- 
sponding equations for the exactly homogeneous perturbations with fc = 0. First note that 
Pi and Qi depend on the metric and the matter variables through their scalar combinations. 
In particular the 4-velocities appear in the form u^uj^ and/or the scalar combinations of 
u^ and its covariant derivatives. In the linear perturbation the terms proportional to vi 
produced from them all come with positive powers of k. This implies that F/ and Egi are 
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written as linear combinations of ^, k'^^, A^/, k'^Z and kVjj = k{vi — vj) with coefficients 
which do not contain negative powers of k. Hence the k ^ hniit of Eq.f p.l?] ) and the 
corresponding k = equation coincide with each other under the condition /c^$ — > 0. 

On the other hand, Eq.( |6.18D , which vanishes identically for k = due to the vector 
origin of F^, gives a non-trivial equation in the /c — *> limit. If we require the condition 



n/-^o, 



it is written as 



hjZj 
H 



Qi 



+ ShiZj + hi A + piTi + c^iPiAgi = ^Z + lim F,i/k. 



fc^O 



(6.20) 



(6.21) 



Since Z is related to A by Eq.( p.5|) , this equation can be regarded as one determining the 
fc — i> limit of Zj, or \iTak^o{vi — B)/k, which has no relation to the quantities describing 
homogeneous perturbations, from a solution for A and Agj to the perturbation equations 
with k = 0(recall the argument on Tj). This is in accordance with the general argument 
in §3. 

In general this gives a set of coupled equations for Zj. However, if the condition 



F,i/k -^ 
is satisfied, it can be explicitly integrated to give 



hiZi = — 



C 



I - I dta^{hjA + piTj + cjpiAgj - Shiq^Z) 
Jt(i 



(6.22) 



(6.23) 



where Ci is an integration constant and to is an initial time. In particular, for the universal 
adiabatic modes given in §4, this gives 



hjZi 



H 



C,^U-K.-\-k. 



(6.24) 



where the subscript denotes the value at t = to- Here note that Eq.( |6.22| ) is a rather 
strong condition in realistic situations where F^i is usually proportional to the relative 
velocities Vjj. 

As discussed in §3, in order for this A; ^ limit solution to correspond to a fc = solu- 
tion, the consistency conditions (|3.3| ) and (|3.4| ) should be satisfied. Under the requirement 
(|6.22| ), these conditions are simply reduced to the asymptotic condition 



kCi{k) -^ 0. 
Further, Z/ is also restricted by the A; — ;► limit equation 
hZ = ^ hjZi 



(6.25) 
Since Z is expressed as 



H 
H 



J2'^i+ [ ^^«^ (-^^ -pT- clpA, 



T.ci- 



3 H 



pA 



H 



pA 



(6.26) 
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from Eq.( |2.25| ) and Eq.( |2.3CI|) , this condition is expressed as 

T^Ci + ll^^pAJ =0. (6.27) 



This is just the equation (|3~5|) at the initial time t = t^. 

For a iV-component perfect fluid system, the perturbation of the exactly homogeneous 
solutions generates A^-independent solutions for Agj and A, which satisfy the fc — > limit of 
Eq. ( |6.17D and the Einstein equations. Each of these solutions in turn determines Zj through 



Eq.( |6.23|) with A^ — 1 independent integration constants satisfying the condition (|6.27|) . By 
this procedure we obtain 2N — 1 independent solutions to the whole perturbation equations 
in the fc — ;> limit. Thus, by adding the universal adiabatic decaying mode corresponding 
to the trivial homogeneous perturbation, we can obtain all the 2N independent k -^ 0-limit 
solutions for this system from the exactly homogeneous solutions. 

§7 Summary and Discussions 

In this paper we have investigated the relation between the fc — ^ limit of solutions to the 
perturbation equations with fc 7^ and the exactly homogeneous perturbations {k = 0) for 
quite general mult i- component systems on a spatially flat Robertson- Walker universe. The 
main result is summarized as follows. 

First for the case in which the fundamental variables describing cosmic matter contain 
a dynamical vector field, for any solution to the perturbation equations with k = there 
exists a solution to the perturbation equations with fc 7^ which coincides with the former 
in the k ^ limit. On the other hand if cosmic matter consists only of scalar fields, 
such a correspondence holds if and only if the exact homogeneous perturbation satisfies 
an additional constraint corresponding to the momentum constraint which is expressed by 
a linear equation on the initial condition for the homogeneous perturbation. This implies 
that for the A/^-component scalar field system one can obtain the k ^ limit of 2N — 1 
independent solutions to the perturbations equations if one knows all the solutions to the 
homogeneous background equations. 

We have also shown that from trivial exactly homogeneous perturbations one can always 
construct two adiabatic solutions to the perturbation equations in the k —>■ limit which 
can be explicitly expressed in terms of integrals of background quantities, irrespective of 
the content and the interactions of cosmic matter. One of them corresponding to the 
vanishing homogeneous perturbation represents a decaying mode for which the Bardeen 
parameter vanishes. The other corresponds to a simple constant scaling of the scale factor 
and gives a growing mode whose Bardeen parameter is a non- vanishing constant. This 
constant scaling is not a gauge transformation, but an extra symmetry of the background 
equations which holds only for spatially homogeneous Robertson- Walker universes. Thus 
this extra symmetry is the hidden reason why exact solutions in the A; — >■ limit have been 
so far found for various systems in spite of the non-trivial structure of their perturbation 
equations. 

For the multi-component scalar field system these universal adiabatic modes and the 
solutions obtained from the non-trivial exact homogeneous perturbations exhaust all the 
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solutions to the perturbation equations in the k ^ hmit. Hence the evolution of pertur- 
bations on superhorizon scales can be determined only from the knowledge on the homoge- 
neous solutions to the background equations in the lowest-order approximation. However, 
if one wants to know the higher-order correction in k, one must solve the perturbation 
equations k ^ iteratively |0. In this procedure one needs the A; — > limit of all the 
solutions to the perturbation equations for the gauge-invariant variable describing the per- 
turbation of the scalar fields. As we have shown, this information can be also obtained 
from the exactly homogeneous solutions to the perturbation equations. Thus the behavior 
of superhorizon perturbations of this system can be determined with any accuracy from 
the knowledge on the homogeneous solutions. 

This argument can be directly extended to the scalar field system coupled with radia- 
tion, for which all the /c — > limit of the solutions to the perturbation equations are directly 
determined from the homogeneous background solutions. This correspondence may pro- 
vide a very powerful method to analyze the evolution of perturbations during reheating in 
the inflationary models because the structure of the equations for the homogeneous scalar 
fields decaying to radiation is much simpler than that for the perturbation equations with 
energy transfer terms. 
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